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We present the conditional quantum dynamics of an electron tunneling between two quantum 
, dots subject to a measurement using a low transparency point contact or tunnel junction. The 

■ double dot system forms a single qubit and the measurement corresponds to a continuous in time 

readout of the occupancy of the quantum dot. We illustrate the difference between conditional 
and unconditional dynamics of the qubit. The conditional dynamics is discussed in two regimes 
depending on the rate of tunneling through the point contact: quantum jumps, in which individual 
electron tunneling current events can be distinguished, and a diffusive dynamics in which individual 
events are ignored, and the time-averaged current is considered as a continuous diffusive variable. 
We include the effect of inefficient measurement and the influence of the relative phase between the 
two tunneling amplitudes of the double dot /point contact system. 
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I. INTRODUCTION 

, In condensed matter physics, measurements are usually made on many identical quantum systems prepared at the 
same time. For example, in nuclear or electron magnetic resonance experiments, an ensemble of systems of nuclei and 
T i electrons are probed to obtain the resonance signals. In this case, the measurement result is an average,response of the 
' O ■ ensemble. On the other hand, various proposed condensed-matter quantum computer architecturealTQ demand the 
5 ' need to readout physical properties of a single electronic qubit, such as charge or spin at a single electron level. This 
is a non-trivial problem since it involves an individual quantum particle measured by a practical detector in a realistic 
environment. It is particularly important to take account of the decoherence introduced by the measurements on the 
qubit as well as to understand how the quantum state of the qubit, conditioned on a particular single realization of 
measurement, evolves in time for the purpose of quantum computing. 

We consider, in this paper, the problem of an electron tunneling between two coupled quantum dots (CQD's), 
a two-state quantum system (qubit), using a low-transparency point contact (PC) or tunnel junction as a detector 
' (environment) continuously mcasurin|S the position of the electron, schematically illustrated in Fig. |l|. This problem 
^f) I has been extensively studied in Refs.Ert3. The case of measurepiettts by a general quantum point contact detector 
' with arbitrary transparency has also been investigated in Refs.ll3"E3. Jp,^add|itjjm,j_^ similar system measured by a 
single electron transistor rather than a PC has been studied in Refs.EllEJa'EiE^Cja. The influence of the detector 
(environment) on the measured system can be determined by the reduced density matrix obtained by tracing out the 
environmental degrees of the freedom in the total, system plus environment, density matrix. The master equation (or 
rate equaiipns for the reduced-density-matrix elements) for the CQD system (qubit) has been derived and analyzed 
in Refs.l3o. However, this (unconditional) master equation describes only the ensemble average property for the 
CQD system. Integrating or tracing out the PC environmental degrees of the freedom is equivalent to completely 
(~| ignoring or averaging over the results of all measurement records (electron current records in this case). In this 
O I sense, the PC is treated as a pure environment rather than a measurement device which can provide information 
O ■ about the change of the state of the qubit. On the other hand, if a measurement is made on the system and the 
^ , results are available, the state or density matrix is a conditional state conditioned on the measurement results. Hence 
the deterministic, unconditional master equation cannot describe the conditional dynamics of the qubit in a single 
realization of continuous measurements which reflects the stochastic nature of an electron tunneling through the PC 
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KorotkovD'ElO has obtained the Langevin rate equations for the CQD system measured by a ideal PC detector. 
These rate equations describe the randorn-evolution of the density matrix thaLboth conditions, and is conditioned 
by, the PC detector output. Recently, wetj presented a quantum trajectory^^t3 measurement analysis of the same 
system. We found that the conditional dynamics of the CQD system can be described by the stochastic Schrodinger 
equation for the conditioned state vector, provided that the information carried away from the CQD system by the 
PC reservoirs can be recovered by the perfect detection of the measurements. We also analyzed the localization rates 
at which the qubit becomes localized in one of the two states when the coupling frequency between the states is 
zero.|-We,siiowed that the localization time discussed there is slightly different from the measurement time defined in 
RefsHMO. The mixing rate at which the two possible states of the qubit become mixed when O 7^ was calculated 
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as well and found in agreement with the result in Ref.Bllil. In this paper, we focus on the qubit dynamics conditioned 
on a particular realization of the actual measured current through the PC device. Especially, we take into account 
the effect of inefficient measurement on the conditional dynamics and illustrate the conditional quantum evolutions 
by numerical simulations. l-l p 

The problem for a "non-ideal" detector was discussed in Refslrlj. There the non-ideality of the detector is modeled 
as two ideal detectors "in parallel" with inaccessible output of the second detector. The information loss is due 
to the interaction with the second detector, treated as a "pure environment" , (which does not affect the detector 
current). As a consequence, the decoherence rate, Ttot, in that case is larger than the decoherence rate for the PC as 
an environment alone, Ttot ~ P^ = 7^ > 0. Hence an extra decoherence term, —jdPab, for example, is added in the 
rate equation pab- However, this approach does not account for the inefficiency in the measurements, which arises 
when the detector sometimes misses detection. In that case, there is still only one PC detector (environment) and 
disregarding all measurement records leads to Ttot = Pd- Furthermore, the detector current is affected and in fact 
reduced by the inefficiency in the measurements. 

In this paper, we analyze the conditional qubit dynamics analytically and numerically. We take into account 
the effect of inefficient measurement of the PC detector on the dynamics of the qubit. The different behavior of 
unconditional and conditional evolution is demonstrated. We pressent the mnditipnal quantum dynamics over the 
full range of behavior, from quantum jumps to quantuift,diffusion.EJ In Refs.BlJ'aEj, the two tunneling amplitudes of 
the CQD-PC model were assumed to be real. In Ref.llJ, the relative phase between them was taken into account. 
Here, we discuss and illustrate furthermore their influence on the qubit dynamics. In Sec. |^ we describe the model 
Hamiltonian and the unconditional master equation. We then obtain in Sec. [II the quantum-jump and quantum- 
diffusive, conditional master equations for the case of inefficient measurements. Sec. |l^ is devoted to the analysis 
for the qubit dynamics. Numerical simulations of the conditional evolution are presented in this section. Finally, a 
short conclusion is given in Sec. In the Appendix, the stationary noise power spectrum of the current fluctuations 
through the PC barrier is calculated in terms of the quantum-jump formalism. 



II. UNCONDITIONAL MASTER EQUATION FOR THE CQD AND PC MODEL 

The appropriate way to approach quantum measurement problems may be to treat the measured systeraL, the 
detector (environment), and the coupling between them together microscopically. Following from Refs.Ell3'El, We 
describe the whole system (see Fig. |^) by the following Hamiltonian: 

— 'HcQD + Ti-pc + 'Hose + Ti-coup (1) 

where 

'HcQD = It- \^ic\ci + UJ2c\c2 + ^{c\c2 + C2Ci)j , (2) 

'Hpc ^^Y^ {^kO,[k'^Lk + LOka^jii^aRk^ + ^ (TfcgaLaflg + Tg^a^j^^aLk^j , (3) 

k k.q 

Hcoup = '^I'^l {xkqO-lk'^Rq + X*qkO-Rq<^Lk'j ■ (4) 

k,q 

H.CQD represents the effective tunneling Hamiltonian for the measured CQD system (mesoscopic qubit). For simplicity, 
we assume strong inner and inter dot Coulomb repulsion, so only one electron can occupy this CQD system. We label 
each dot with an index 1, 2 (see Fig. |l|) and let Ci (cl) and huji represent the electron annihilation (creation) operator 
and energy for a single electron state in each dot respectively. The coupling between these two dots is given by hfl. 
The tunneling Hamiltonian for the PC detector is represented by Tipc where ULk, a.Rk and hujj^, hcu^ are respectively 
the fermion (electron) field annihilation operators and energies for the left and right reservoir states at wave number k. 
One should not be confused by the electron in the CQD with the electrons in the PC reservoirs. The tunneling matrix 
element between states k and q in left and right reservoir respectively is given by Tkq. Eq. (Q), Ticoup, describes the 
interaction between the detector and the measured system, depending on which dot is occupied. When the electron 
in the CQD system is near the PC (i.e., dot 1 is occupied), there is a change in the PC tunneling barrier. This barrier 
change results in a change of the effective tunneling amplitude from Tkq Tkq + Xkq- As a consequence, the current 
through the PC is also modified. This changed current can be detected, and thus a measurement of the location of 
the electron in the CQD system is effected. 

The (jjpconditional) master equation of the reduced density matrix ijax the CQD system (qubit^ has been obtained 
in RefsH'til. Under similar assumptions and approximations as in Ref.EJ, the zero-temperaturej^ Markovian master 
equation of the qubit can be written as: 
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p{t) = -'j^[ncQD,p{t)] + V[r + Xn^]p{t) (5a) 
= Cp{t), (5b) 

where ni = c|ci is the occupation number operator for dot 1 and the parameters T and X are given by 

irp ^D = 27re\Too\^ gL9RV/h, (6a) 
\T + ^D' = 27re|roo + Xm?gL9RVin, (6b) 

where D and I?' are the average electron tunneUng rates through the PC barrier without and with the presence of 
the electron in dot 1 respectively. In Eq. (H), eV — p,L ~ P-bAs the external bias applied across the PC (where e is 
the electronic charge, and pL and pn stand for the chemical potentials in the left and right reservoirs respectively), 
Too and Xoo are energy-independent tunneling amplitudes near the average chemical potential, and gL r^fMr#R are 
the energy-independent density of states for the left and right reservoirs. In Eq. (|5a|), the superoperatorE3'E3E3 V is 
defined as: 

V[B]p^J[B]p-A[B]p, (7) 

where 

J[B]p = BpB\ (8) 
A[B]p:^ {B^Bp + pB''B)/2. (9) 

Finally, Eq. ( |5b| ) defines the Liouvillian operator C. The form of the master equation (|5^), defined through the 
superoperator 'D[B]p{t), preserves the positivity of the density matrix operator p{t). Such a Markovian master 
equation is called a LindbladEO form. P 

Evaluating the density matrix operator in the same basis as in Ref.El, we obtain 

Paa{t) = i^[pab{t) - Pba{t)] , (10a) 
Pabit) = iSpabit) + in[paait) - pbb{t)] - (| A"] V2)pah(t) + i lm{T* X) Pabit) , (10b) 

where h£ — h{uj2 — wi) is the energy mismatch between the two dots, Pij{t) — {i\p{t)\j), and Paa{t) and pbb{t) are the 
probabilities of finding the electron in dot 1 and dot 2 respectively. The corresponding logical qubit states are then 
I a) and \b). The rate equations for the other two density matrix elements can be easily obtained from the relations: 
Pbb{t) = 1 — paa{t) and pba{i) = Pabi^) ■ Compared to an isolated CQD system, the presence of the PC detector 
(environment) introduces two effects to the CQD system. First, the imaginary part of {T*X) (the last term in Eq. 
( |lOb| )) causes an effective shift in the energy mismatch between the two dots. Second, it generates a decoherence 
(dephasing) rate 

rd = \x\y2 (11) 

for the off-diagonal density matrix elements. The relative phase between the two complex tunneling amplitudes may 
produce additional effects on conditional dynamics of the CQD system as well. This will be shown later when we 
discuss conditional dynamics. Physically, the presence of the electron in dot 1 (state |a)) raises the effective tunneling 
barrier of the PC due to electrostatic repulsion. As a consequence, the effective tunneling amplitude becomes lower, 
i.e., D' = \T + X\'^ < D = \T\'^. This sets a condition on the relative phase 9 between A" and T: cos 6* < ~'\X\/{2\T\). 



III. CONDITIONAL MASTER EQUATION 

Equation (^ describes the time evolution of reduced density matrix when all the measurement results are ignored, 
or averaged over. To make contact with a single realization of the measurement records and study the stochastic 
evolution of the quantum state, conditioned on a particular measurement realization, we obtain in this section the 
conditional master equation. 

The measurable quantities, such as accumulated number of electrons tunneling through the PC barrier, are stochas- 
tic. On average of course the same current flows in both reservoirs. However, the current is actually made up of 
contributions from random pulses in each reservoir, separated in time by the times at which the electrons tunnel 
through the PC, assuming the response of the reservoirs is fast. In this section, we first treat the electron timneling 
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current consisting of a sequence of random 5 functirm. pulses. In other words, the measured current is regarded as a 
series of point processes (a quantum-jump model)E3H. The case of quantum diffusion will be analyzed later. 

For almost-all infinitesimal time intervals, the measurement result is null (no detection of an electron tunneling 
through the PC barrier). The system in this case changes infinitesimally, but not unitarily. The nonunitary component 
reflects the changing probabilities for future events conditioned on past null events. At randomly determined times 
(conditionally Poisson distributed), there is a detection result. When this occurs, the system undergoes a finite 
evolution, called a quantum jump. 

Formally, we can write the current through the PC as 



iit) ^edN{t)/dt, 



(12) 



where dN(t) is a classical point process which represents the number (either zero or one) of Uumeling events seen in 
an infinitesimal time dt. We can think of dN{t) as the incremftst in the accumulated numberO'EilEil of electrons N{t) 
in the drain in time dt. The point process is formally definedlil by the conditions on the classical random variable 
dNcit): 



[dN,{t)]^ ^dN,{t), 
E[dN,{t)]=CTT[pi,{t- 



dt)] =([D+{D' - D){ni)c{t)]dt. 



where {ni)c{t) = Tr[ni/9c(t)], E\Y] denotes an ensemble average of a classical stochastic process F, and 

pic{t + dt) = J\r + Xn^\prXt)dt 



(13a) 
(13b) 



(14) 



is the unnormalized density matrixtil given the result of an electron tunneling through the PC barrier at the end of 
the time interval -I- dt). We explicitly use the subscript c to indicate that the quantity to which it is attached is 
conditioned on previous measurement results, the occurrences (detection records) of the electrons tunneling through 
the PC barrier in the past. Note that the density matrix Pc(t) is not the solutio n of the miconditional reduced master 
equation, Eq. (|5a|). It is actually conditioned by dN^it!) for t' < t. Equat ion ( 13a ) simply states that dNc{t) equals 
either zero or one, which is why it is called a point process. Equation (13b) indicates that the ensemble average 
of dNc{t) equals the probability (quantum average) of detecting electrons tunneling through the PC barrier in time 
dt. The inefficiency in the measurements, which arises when the detector sometimes misses detection, is taken into 
account here. The factor C ^ 1 represents the fraction of detections which are actually registered by the PC dete ctor . 
The value C = 1 then corresponds to a perfect detector or efficient measurement. By using Eq. (p^), equation (p^lj 
with C = 1 states that the average current is eD when dot 1 is empty, and is eD' when dot 1 is occupied. In RefH 
the case of inefficient measurements is discussed in terms of insufficiently small readout period. In other words, the 
bandwidth of the measurement device is not large enough to resolve and record every electron tunneling through the 
PC barrier. I— ■ 

By following the similar derivation as in Ref.til, the stochastic master equation of the density matrix operator, 
conditioned on the observed event in the efficient measurement in time dt can be obtained: 



dpc{t) = dN^t) 



J[T + Xni\ 

Vlc{t) 



- 1 



Pc{t) 



+ dt \ -A[T + Xni]p,{t) + (1 - Q)J[T + Xni\p,{t) + C'Pic{t)pc{t) + jXHcqd, Pc{t)] 



(15) 



where 



Vic{t) = D + {D' - D){ni)c{t). 



(16) 



Aver aging this equation over the observed stochastic process, by setting E[dNc{t)] equal to its expected value Eq. 
( 13t ) , gives the unconditional, deterministic master equation (|5a|) . |— ■ 

Similarly, the extension to the case of quantum diffusion can be carried out as in Ref.EJ. In this case, the average 
electron tunneling rate is very large compared to the extra change of the tunneling rate due to the presence of the 
electron in the dot closer to the PC. In addition, individual electrons tunneling through the PC are ignored and time 
averaging of the currents is performed. This allows electron counts, or accumulated^electron number, to be considered 
as a continuous diffusive variable satisfying a Gaussian white noise distributiorO'L3 in time 5t: 



6N{t)^{(:\T\'[l + 2ecose{ni),{t)] + ^Q \r\i{t)}5t, 



(17) 
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where e = (|A'|/|T|) ^ 1, 6 is the relative phase between X and T, and ^{t) is a Gaussian white noise characterized 
by 



Here E denotes an ensemble average and S{t — t') is a delta function. In stochastic calculusE3t3, (,{t)dt = dW{t) is 
known as the infinitesimal Wiener increment with the property E[{(,(t)dty] = E[dW{t)'^] = dt. In obtaining Eq. dlj), 
we have assumed that 2|T||A'| cos9 ^ \X\'^. Hence, for the quantum-diffusive equations obtained later, we should 
regard, to the order of magnitude, that | cos0| ^ 0{\) ^ e = (jA'l/ITI) and | sin 6*1 ^ 0{e) ^ 1. The quantum-diffusive 
conditional master equation can be found as: 

Pc{t) =-^[HcQD,Pc{t)]+V[T + Xni]p,{t) 



<it)^ [r*Xnip,{t) + X*Tp,{t)ni - 2 Re(r*A')(ni),(i)p,(i)] . (19) 



In arriving at Eq. (pD|), we have used the stochastic Ito calculust2rE3 for the definition of derivative as p{t) = 
limdt^o[p(^ + dt) — p(t)]/dt. This is in contrast to the defiaition, p{t) = limdt^o[p{t + dt/2) — p{t — dt/2)]/dt, used 
in another stochastic calculus, the Stratonovich calculua3c3. It is easy to see that the ensemble average evolution 
of Eq. ( [l9| ) reproduces the unconditional master equation (|5^) by simply eliminating the white noise term using Eq. 

If and only if detections are perfect (efficient measurement), i.e. C — 1; the stochastic master equations for the 
conditioned density matrix operators, ( |l5|) and (|is|), equivalent to the stochastic Schrodinger equations (Eqs. (35) 
and (41) of Ref.ll3, respectively) for the conditioned states. That is because in order for the system to be continuously 
described by a state vector (rather than a general density matrix), it is necessary (and sufficient) to have maximal 
knowledge of its change of state. This requires perfect detection or efficient measurement, which recovers and contains 
all the information lost from the system to the reservoirs. If the detection is not perfect and some information about 
the system is unrecoverable, the evolution of the system can no longer be described by a pure state vector. For the 
extreme case of zero efficiency detection, the information (measurement results at the detector) carried away from 
the system to the reservoirs is (are) completely ignored, so that the stochastic master equations (|l^) and ([l^) after 
being averaged over all possible measurement records reduces to the unconditional, deterministic master equation 
(^), leading to decoherence for the system. 

To make the quantum-diffusive, conditional stochastic master equation (|l^ more transparent, we evaluate Eq. ( [l9| ) 
in the same basis as for Eq. ( |lo|) and obtain: 



Paait) = in[pab{t) - pfc,(t)] + ^8 C Trf COS 9 Paa{t) Pbbit)^ , (20a) 
Pabit) = l{£ + \T\\X\ Sine)pabit) + tn[paait) - Pbbit)] - Erf p,b(i) 

+ ^2 CT d{cos 9 [pbbit) ~ Paait)] +tsine} pab{t)m , (20b) 

where we have set | A"! = \/2rrf. Again, the ensemble average of Eq. ( pO| ) by eliminating the white noise terms reduces 
to Eq. (^. It is also easy to verify that for zero efficiency C = 0, the conditional equations ([l5|), (pj|), and (|2^), 
reduce to the corresponding unconditional ones, ( pa| ) and ( p^ ) respectively. That is, the effect of averaging over all 
possible measurement records is equivalent to the effect of completely ignoring the detection records or the effect of 
no detection results being available. 

IV. CONDITIONAL DYNAMICS UNDER CONTINUOUS MEASUREMENTS 

In this section, we analyze the qubit dynamics in detail and present the numerical simulations for the time evolution. 
We represent the qubit density matrix elements in terms of Bloch sphere variables since some physical insights into 
the dynamics of the qubit can sometimes be more easily visualized in this representation. Denoting 

I = clc2+c\ci, (21) 

ax = eld + c\c2, (22) 

ay = —ic\ci + ic\c2, (23) 

o-z = - c|ci, (24) 
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we write the density matrix operator for the CQD system (qubit) in terms of the Bloch sphere vector {x,y, z) as: 



p{t) = [I + x{t)a^ + y{t)ay + z{t)a,]/2 



1 + z{t) x{t) - iy{t) 
x{t)+iy{t) l-z{t) 



(25a) 
(25b) 



It is easy to see that Trp(t) = 1, / is a unit operator, Ui defined above satisfies the properties of Pauh matrices, 
and the averages of the operators Ux, Uy, are a;(t), j/(t), z(t) respectively. In this representation, the variable 
z(t) represents the population difference between the two dots. Especially, z(t) — 1 and z(t) = —I indicate that the 
electron is localized in dot 2 and dot 1 respectively. The value z(t) = corresponds to an equal probability for the 
electron to be in each dot. Generally the product of the off-diagonal elements of p{t) is smaller than the product of 
the diagonal elements, leading to the relation x'^{t) + y^{t) + z^{t) < 1. When p{t) is represented by a pure state, the 
equal sign holds. In this case, the system state can be characterized by a point {x,y, z) on the Bloch unit sphere. 

The master equations (|5a|), ( pj| ) and (p^), can be writt en as a set of coupled stochastic differential equations in terms 
of the Bloch sphere variables. By substituting Eq. (25a) into Eq. (|5a|), and collecting and equating the coefficients in 
front of ax, fJy, respectively, the unconditional master equation under the assumption of real tunneling amplitudes 
is equivalent to the following equations: 



dxit) 
dt 

dyjt) 

dt 
dz{t) 
dt 



^ -{£ + \T\\X\sine)y{t) - Vdx{t), 

= {E + \T\\X\ sin9)x{t) - 2^z{t) - Tayit), 

= 2ny{t). 



Similarly for the quantum-diffusive, conditional master equation (|T^), we obtain 

dxc{t) 



dt 
dycjt) 

dt 
dzcjt) 

dt 



-{£ + \T\ \X\ sin 6) y,{t) - Td x,{t) + ^2CTd[- sin 9 ydt) + cos d Zc{t)xc{tMt), 



{£ + \T\ \X\ sin 6) xc{t) ~ 217 z,{t) - ydt) + ^2CTd[sm 9 x^t) + cos 9 zdt)ycitMt), 
2nydt) - V2CT~d cos 9 [l - z^Jt)] ^t). 



(26a) 
(26b) 
(26c) 

(27a) 
(27b) 
(27c) 



Again the c-subscript is to emphasize that these variables refer to the conditional state. It is trivial to see that Eq. 
( p7| ) averaged over the white noise reduces to Eq. (p^), provided that E[xc{t)] — xit) as well as similar replacements 
are performed fo r vA f) and Zc{t). The analogous calculation can be carried out for the quantum-jump, conditional 
master equation (|l^). We obtain 



dxdt) =dt[-[£ + {l^ 0\'T\\X\ sin9] y^t) - (1 - C)Td x^t) 



C {D' - D) 



zdt)Xc{t) 



-dNdt) 



2\T\\X\ sm9ydt) + [2rrf - jP' - D)zdt)]xdt) 
2D + {D' - D)[l ~ zc{t)] 



dydt) = dt 



[£ + {1- C)\T\\X\ sin9] xdt) - (1 - O^dVdi) - inz,{t) - ^K dt)ydt) 



-dNdt) 



-2\T\\X\ sm9xdt) + [2rd - [D' - D)zdt)]ydt) 



dzdt) = dt 2^1 ydt) 



2D + {D' - D)[l - zdt)] 
C{D'-D) 



[l-^cW])-rf^c(t)( 



{D'-D)[l-zm 
2D + {D' - D)[l ~ zdt)] 



(28a) 

(28b) 
(28c) 



As expected, by using Eq. ( 13b ), the ensemble average of Eq. (E8h also reduces to the unconditional equation (|2q). 
One can also observe that for ( ~ 0, the conditional equations, (p8|) and (p7|), reduce to the unconditional equation 
(Efl)_as well. Next, we analyze the conditional qubit dynamics. Part of the results in Sec. IV A have been reported in 
Rcf.y. 
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A. From quantum jumps to quantum diffusion 



Fig. ^a) shows the unconditional (ensemble average) time evolution of the population difference z{t) with the initial 
qubit state being in state |a), i.e., dot 1 is occupied . The unconditional population difference z(t), rises from — 1, 
undergoing some oscillations, and then tends towards zero, a steady (maximally mixed) state. On the other hand, 
the conditional time evolution, conditioned on one possible individual realization of the sequence of measurement 
results, behaves quite differently. We consider first the situation, where D' = \T + X]"^ = 0, discussed in Ref.cl. In 
this case, due to the electrostatic repulsion generated by the electron, the PC is blocked (no electron is transmitted) 
when dot 1 is occupied. As a consequence, whenever there is a detection of an electron tunneling through the PC 
barrier, the qubit state is collapsed into state \b), i.e., dot 2 is occupied. The quantum-jump conditional evolution 
shown in Fig. ^(b) (using the same parameters and initial condition as in Fig. ^(a)) is rather obviously different from 
the unconditional one in Fig. ^(a). The conditional time evolution is not smooth, but exhibits jumps, and it does not 
tend towards a steady state. One can see that initially the system starts to undergo an oscillation. As the population 
difference Zc{t) changes in time, the probability for an electron tunneling through the PC barrier increases. This 
oscillation is then interrupted by the detection of an electron tunneling through the PC barrier, which bring Zc{t) 
to the value 1, i.e., the qubit state is collapsed into state \b). Then the whole process starts again. The randomly 
distributed moments of detections, dNc(t), corresponding to the quantum jumps in Fig. ||(b) is illustrated in Fig. 
^(c). Although little similarity can be observed between the time evolution in Fig. ^(a) and (b), averaging over many 
individual realizations shown in Fig. |^(b) leads to a closer and closer approximation of the ensemble average in Fig. 
1(a). 

Nexipjire illustrat e h ow the transition from the quantum-jump picture to the quantum-diffusive picture takes place. 
In ReiM and Sec. Ill, we have seen that the quantum-diffusive equations can be obtained from the quantum-jump 
description under the assumption of |T| 3> \X\. In Fig. ^(a)-(d) we plot conditional, quantum-jump evolution of 
Zc{t) and the corresponding moments of detections dNc{t), with different (|T|/|A'|) ratios. Each jump (discontinuity) 
in the Zc{t) curves corresponds to the detection of an electron through the PC barrier. One can clearly observes 
that with increasing (|T|/|A'|) ratio, the number of jumps increases. The amplitudes of the jumps of Zc{t), however, 
decreases from D' = with the certainty of the qubit being in state \b) to the case of {D — D') « {D -\- D') with 
a smaller probability of finding the qubit in state Nevertheless, the population difference Zc{t) always jumps up 
since D = |Tp > D' = \T + X\'^. In other words, whenever there is a detection of an electron passing through PC, dot 
2 is more likely occupied than dot 1. The case for quantum diffusion using Eq. ( |27| ) is plotted in Fig. ||(e). In this case, 
very small jumps occur very frequently. We can see that the behavior of Zc{t) for |T| = 5\X\ in the quantum-jump 
case shown in Fig. |3|(d) is already very close to that of quantum diffusion shown in Fig. ^(c). To minimize the number 
of controllable variables, the same randomness is applied to produce the quantum-jump, conditional evolutions in Fig. 
^(a)-(d). This, however, does not mean that they would have had the same detection output, dNc{t). The number 
of tunneling events in time dt, dN^it), does not depend on the randomness alone. It also depends on |T|, and 
9, and has to satisfy Eq. ( |l3bD in a self-consistent manner. In fact, it both conditions and is conditioned by the 
conditional qubit density matrix. Note that the unconditional evolution does not depend on the parameter |T| when 
6 = TT (see Eq. (26)). This implies that depending on the actual measured detection events, different measurement 
schemes (measurement devices with different tunneling barriers or different values of |T| when 9 = ir) give different 
conditional quantum evolutions. But they would have the same ensemble average property if other parameters and 
the initial condition are the same. Hence, averaging over all possible realizations, for each measurement scheme in 
Fig. H, will lead to the same ensemble average behavior shown in Fig. ||(a). 



B. Quantum Zeno effect 



The quantum Zeno effeclxan be naturally described by the conditional dynamics. The case for quantum diffusion 
has been discussed in Ref.Llu. Here, for completeness, we discuss the quantum-jump case. The quantum Zeno effect 
states that repeated observations of the system slow down transitions between quantum states due to the collapse of 
the wave function into the observed state. Alternatively, the interaction with one measurement apparatus destroys 
the quantum coherence (oscillations) between \a) and |&) at a rate that is much faster than the tunneling rate fl. 
For fixed fi, |T| and 9, by increasing the interaction with the PC detector \X\ = y/2Td, we increase the number 
and amplitude of jumps and hence the probability of the wave function being collapsed to the localized state. The 
time evolutions of the population difference Zc{t) for different ratios of (F^/ri) are shown in Fig. ^. Here, the initial 
qubit state is |a), and other parameters are: ( = 1, £ ~ 0, 9 = tt, |Tp = lOfl. We can observe that the period of 
coherent oscillations between the two qubit states increases with increasing (r^/ri), while the time of a transition 
(switching time) decreases. In the limit of vanishing f2, a transition from one qubit state to the other state takes a 
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time (switching time) of order of locaUzation time0, I/tj'™'' = {D + D')/\rd{y/ D + \/ D')"^]. In the parameter regime 
of Fig. ^c) (Fcf/fi = 8), this time is stiU much smaller than the average time between state-changing transitions 
(period of oscillations) due to il, i.e., the mixing timelia, l/7mix = rrf/(4ri^). Hence, we can already see from Fig. 
^(c) for Fd/f^ = 8 that very frequent repeated measurements would tend to localize the system. 

The ensemble average behavior of z{t) is also shown in dashed line in Fig. ^ E — and initially the electron is 
in dot 1, from the solution of Eq. (|lo| ) , the probability Paa{t) = [1 ~ ^(^)]/2 can be written as 



cosh ( I H ^ sinh ( —^t 



(29) 



where fir = \/F^ — (417)2. ^j^g Appendix, the stationary noise power spectrum jaf the current fluctuations through 
the PC barrier is calculated for the case of £ = and the result can be written asO: 



Situ) = 5o + 



(cj2 _ 4f]2)2 _^ 



(30) 



where So — 2eioo — C {D' + D) represents the shot noise, ioo = ^({D' + D)/2 is the steady-state current and 
Ai = e({D — D') represents the difference between the two average currents. For F^ < 417, paa{t) shows the damped 
oscillatory behavior in the immediate time regime (see dashed line in Fig. ^(a) and (b)). In this case, the spectrum 
has a double peak structure, indicating that coherent tunneling is taking place between the two qubit states. This 
is illustrated in Fig. ||(a) and (b). When Vd > 417, paa{t) does not oscillate but decays in time purely exponentially, 
saturating at the probability 1/2. (see dashed line in Fig. ^(c)) This corresponds to a classical, incoherent behavior. 
In this case, only a single peak, centering at = 0, appears in the noise spectrum, as illustrated in Fig. |^(c). The 
evolution of Zc(t) in Fig. ^c), is one of the possible conditional evolutions in this parameter regime (Fd/li — 8). 
In this parameter regime F^ > 417, the conditional evolution Zc{t) behaves very close to a probabilistic jumping or 
random telegraph process. After ensemble averaging over all possible realizations of such conditional evolutions, one 
would then obtain the classical, incoherent behavior. 



C. Relative phase of the tunneling amplitudes 

The relative phase between the two complex tunneling amplitudes produces effects on both conditional and uncon- 
ditional dynamics of the qubit. In the following, we consider the case that C, — 1 and E — Q. From Eq. (|2^), after each 
jump the imaginary part of the product [T* X) seems to cause an additional rotation around the z-axis in the Bloch 
sphere, but does not directly change the population probability Zc{t) of the qubit. However, the actual conditional 
evolution of the Bloch sphere variables is complicated. It is stochastic and nonlinear, and depends on the relative 
phase of the tunneling amplitudes in a nontrivial way. Nevertheless, after ensemble average, the imaginary part of 
{T*X) generates an effective shift in the energy mismatch of the qubit states (see Eq. (p^). 

There are situations in which the effect of the relative phase of the tunneling amplitudes can be easily seen. For 
C = 1 and f = 0, if the tunneling amplitudes are real, i.e., 9 = 7:, and the initial condition a;c(0) = 0, then the time 
evolution of Xc{t), from Eq. (p8|), does not change and remains at the value at all times. But if 6* ^ tt or smO ^ 0, 
the conditional evolution of Xc{t) behaves rather differently. It changes after the first detection (quantum jump) takes 
place. Fig. ^ shows the evolutions of the Bloch variables Xc(i), yc{i): Zc{t) with the same initial condition (the qubit 
being in |a)) and parameters but different relative phases: = tt for (a)-(c), and 9 = cos^^{\X\/\T\) for (d)-(f). We 
can clearly see quite different behaviors of Xc{t) in these two cases. The asymmetry of the electron population in Zc{t), 
due to effectively generated energy mismatch in the second case in Fig. ^(f), can be roughly observed. The effect of 
the relative phase is small in the case of quantum diffusion. As noted in Sec. |l|, in order for the quantum-diffusive 
equations to be valid, we should regard, to the order of magnitude, that |cos6'| ^ 0(1) and |sin0| ^ 0{e). This 
implies that in this case 9 n. Hence the effect of the relative phase is small and_th£_jConditional dynamics does not 
deviate much from the case that the tunneling amplitudes are assumed to be real3'LfE3. 



D. Inefficient measurement and non-ideality 

We have shownQ that for C = 1, the conditional time evolution of the qubit can be described by a ket state vector 
satisfying the stochastic Schrodinger equation. It is then obvious that perfect detection or efficient measurement 
preserves state purity for a pure initial state. However, the inefficiency and non-ideality of the detector spoils this 
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picture. The decrease in our knowledge of the qubit state leads to partial decoherence for the qubit state. We next 
find the partial decoherece rate introduced in this way. . .■ . 

The stochastic differential equations in the form of Ito calculusc§E^ have the advantage that it is easy to see that 
the ensemble average of the conditional equations over the random process ^(i) leads to the unconditional equations. 



However, it is not a natural physical choice. For example, for ( — I, the term —TdPab{t) in Eq. (20b) does not really 
cause decoherence of the conditional qubit density matrix. It simply compensates the noise term due to the definition 
of derivative in Ito calculus. Hence, in this case the conditional evolution of pab{t) does not really de creas e in time 



exponentially. To find the psiiial decoherence rate generated by inefficiency C < 1, we transform Eq. (20t) into the 
form of Stratonovich calculusHEJ. We then obtain for 9 — tt: 



Pab{t) = i£pab{t) + iil[paa{t) - Pbb{t)] - [pbbif) " Paa{t)] [i{t) ~ Iq] Pab{t) - (1 - O^d Pab{t), (31) 




where i{t) — io = e\T\{( ^/2^d[l — 2paait)] + v^^(^)|. Here we have used the following relations: the conditional 
current i{t) = e5N{t)/5t with 5N{t) given by Eq. (|7| and the average current io = eC(-D + D')/2, where D = |Tp 
and D' = |Tp — 2|T||A'| in the quantum-diffusive limit. In this form, Eq. ( |3l| ) elegantly shows how the qubit density 
matrix is conditioned on the measured current. We find that the last term in Eq. ( |3l| ) is responsible for decoherence. 
In other words, the partial decoherence rate for an individual realization of inefficient measurements is (1 — Cl^d- 
For a perfect detector C = 1; t^^is decoherence rate vanishes and the conditional Pab{t), as expected, does not decay 
exponentially in time. Similar conclusion could be drawn from Eq. ( p8| ) for the quantum-jump case. For 9 = tt, the 
off-diagonal variables Xc{t) and ydt) seem to decrease in time with the rate (1 — ()Td- 

In Bloch sphere variable representation, we can use the quantity Pc{t) — x'^{t) + y1{t) -f z1{t) as a measure of the 
purity of the qubit state, or equivalently as a measure of how much information the conditional measurement record 
gives about the qubit state. If the conditional state of the qubit is a pure state then Pc{t) = 1; if it is a maximally 
incoherent mixed state then Pc{t) = 0. We plot in Fig. ^the quantum-jump, conditional evolution of the purity Pc{t) 
for different inefficiencies, C = li 0.6, 0.2 (in solid linej7 and (in dotted line). Fig. ^(a) is for a initial qubit state 
being in a pure state |a), while Fig. ^(b) is for a maximally mixed initial state. We can see from Fig. 0(a) that the 
purity Pc{t) = 1 at all times for C = 1: while it hardly or not at all reaches 1 for almost all time for ^ < 1. This 
means that partial information about the changes of the qubit state is lost irretrievably in inefficient measurements. 
In addition, roughly speaking, the overall behavior of Pc{t) decreases with decreasing C,. This indicates that after 
being averaged over a long period of time, {Pc{t))t would also decrease with decreasing C. For = 0, the evolution 
of P(t) becomes smooth and tends toward the value zero (the maximally mixed steady state). For a non-pure initial 
state (see Fig. @(b)), the qubit state is eventually collapsed towards a pure state and then remains in a pure state for 
C = 1. But the complete purification of the qubit state cannot be achieved for C < 1- As in Fig. ||(a)-(d), the same 
randomness has been applied to generate the quantum-jump, conditional evolution in Fig. 0(a) and (b). Note that 
the only difference between evolution in Fig. 0(a) and the corresponding one in Fig. 0(b) is the different initial states. 
So when the qubit density matrix in Fig. 0(b) gradually evolves into the same state as in Fig. 0(a), the corresponding 
Pc{t) in Fig. 0(b) would then follow the same evolution as in Fig. 0(a). This behavior can be observed in Fig. 0. The 
purity-preserving conditional evolution for a pjuije iaitial state, and gradual purification for a non-pure initial state for 

atum-diffusive limit. 

|12 



another ideal detector "in parallel" to the original 
lO, an extra-|te|iaM, —'jdPabit), to Eq. ( ^l|) to account 



an ideal detector have been discussed in RefsE'S-El in the 
The non-ideality of the PC detector is modeled in Refs.l 
one but with inaccessible output. We can add, as in Refs.L" 

for the "non- ideality" of the detector. The ideal factor 77 introduced thereCTEIIO can be modified to take account of 
inefficient measurement discussed here. We find 

^ = 1-^ = ^.^, (32) 

i-tot i- d + Jd 

where F = (1 — QTd + jd and Ttot = ^d + Id- For 7^ — 0, we have rj = In RefEl, inefficient measurement is 
discussed in terms of insufficiently small readout period. As a result, the information about the tunneling times of 
the electrons passing through the PC barrier is partially lost. 



V. CONCLUSION 



We have obtained the quantum-jump and quantum-diffusive, conditional master equations for the case of inefficient 
measurements. These conditional master equations describe the random evolution of the measured qubit density 
matrix, which both conditions and is conditioned on, a particular realization of the measured current. We have 
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analyzed the conditional qubit dynamics in detail and illustrated the conditional evolution by numerical simulations. 
Specifically, the conditional qubit dynamics evolving from quantum jumps to quantum diffusion has been presented. 
Furthermore, we have described the quantum Zeno effect in terms of the quantum-jump, conditional dynamics. We 
have also discussed the effect of inefficient measurement and the influence of relative phase between the two tunneling 
amplitudes on the qubit dynamics. 
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APPENDIX A: CALCULATION OF THE NOISE POWER SPECTRUM OF THE CURRENT 

FLUCTUATIONS 

In this Appendix, we calculate the stationary noise power spectrum of the current fluctuations through the PC when 
there is the possibility of coherent tunneling between the two qubit states. Usually one can calculate this noise power 
spectrum using the unconditional, determistic master equation approach, which gives only the average characteristics. 
We, however, calculate it through the stochastic formalism presented here. The fluctuations in the observed current, 
i{t), are quantified by the two-time correlation function: 

G(r) = E[i{t + T)i{t)] - E[i{t + r)] E[i{t)] (Al) 

The noise power spectrum of the current is then given by 

/oo 
dT G{t) e-"^'' . (A2) 
-oo 

The ensemble expectation values of the two-time correlation function for the current in the case of quantum diffusion 
has been calculated in Ref.Ll. Here we will present the quantum-juiBB case. The current in this case is given by Eq. 
(O) . We will follow closely the calculation in the Appendix of Ref .cll to calculate the two-time correlation function, 
E\dNc{t + T)dN{t)]. First we consider the case when t ^ dt > 0, where dt is the minimum time step considered. 
Since dN{t) is a classical point process, it is either zero or one. As a result, E[dNc{t + T)dN{t)] is non-vanishing only 
if there is electron-tunneling event inside each of these two infinitesimal time intervals, [t, t + dt] and [t + T,t + t + dt]. 
Hence, we can write 

E[dNc{t + T)dN{t)] = Proh[dNit) = 1] E[dN4t + T)UNit)=i], (A3) 



where the subscript to the vertical line is the condition for which the subscript on dNdt+r) exists. From Eqs. ( |13b| ) and 



(H), we have Vmh[dN{t) = 1] = C Tr[/5i(i dt]! and E^ 
Using the fact that E[pc{t)] — p{t) and Eqs. ( |5b| ) and 



dN,{t + T)\dN(t)=i]=C^^{J[r + Xni]E[pic{t + T)UN(t)=i]}- 
14f), we can write 



E[pic{t + T)UN[t)=i] = e^(^-'^*)pi(i + dt)/TT[p^{t + dt)] 

= C e'^^^-"''\j[T + Xni]p{t)dt}/'Tv[pi{t + dt)]. (A4) 

Hence, to leading order in dt, we obtain for r > 0: 

E[dNc{t + T)dN{t)] = C^dt^rpj. [ ^ Xni]e'^^{J[T + Xni]p{t)}] . (A5) 

For r = 0, we have, from Eq. (|l3|), that: 

E[dN{t)dN{t)] ^ E[dN{t)] ^C,[D+{D' - D){ni){t)]dt. (A6) 

For short times, this term dominates and we may regard dN{t)/dt as (5-correlated noise for a suitably defined S 
function. Thus the current-current two-time correlation function for r > can be written as: 
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= C + {D' - D)Tv[nip{t)]} 5{t) + C' Tr [j[T + Xni]e^^ {J[T + Xni]p{t)}] . (A7) 

In this form, we have related the ensemble averages of classical random variable to the quantum averages with respect 
to the qubit density matrix. The case t < is covered by the fact that the current-current two-time correlation 
function or G(t) is symmetric in r, i.e., G(r) = G(— r). 



Next we calculate steady-state G(r) and S{uj). We can simplify Eq. (A7) using the following identities for an 
arbitrary operator B: Tr [J[ni]B] = Tr[niB], Tr[e^''B] = Tt[B], and Tr[Be^oo] = Ti[Bp^], where the oo subscript 
indicates that the system is at the steady state and the steady-state density matrix p^c is a maximally mixed state. 
Hence we obtain the steady-state G(r) for r > as: 

G{T)^ei^S{T)+e^C^D'-Df{TT[n,e^^[n,p^]-TT[n,p^]^}, (A8) 



where the steady-state average current ioc — e-CiD + D')/2. The first term in Eq. ( |A8[ ) represents the shot noise 
component. It is easy to evaluate Eq. ( |A8| ) analytically for £ — case. The case for the asymmetric qubit, £ 0, 
can be calculated numerically. Evaluating Eq. ( [A^ ) for £ = 0, wc find 

G{T)=eiooS{T) + ^—^[^^ , A9 

4 V M+ - M- / 



where p± = —{Td/2) ± (rd/2)2 — 4f72, and we have represented Ai — e( {D — D') as the difference between the two 



average currents. After Fourier transform following from Eq. (A2), the power spectrum of the noise is then obtained 
as the expression of Eq. (^). Note that from Eq. (|30|), the noise spectrum at w = 20 for 9 = i.e., real tunneling 
amplitudes, can be written as: 

S{m-So_ WdWW .^10) 
So {D + D') ^^^^> 



where 5*0 = 2eioo = C,{D' + D) represents the shot noise. In obtaining Eq. (AlO), we have used the relation 
Td = {VD — •\/D')^/2 for the case of real tunneling amplitudes. In the quantum-diffusive limit |T| ^ \X\ or 
{D + D') > (D - D'), this ratio [S{2n) - So]/So 4C, independent^ of the values |Qf fl and r^. These results for 
C = 1 and in the limit of quantum diffusion are consistent with those derived in Ref.EI using both the unconditional 
master equation approach and conditional stochastic formalism with white noise current fiuctuations for an ideal 
detector. 
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CQD 




dot 2 



FIG. 1. Schematic representation of an electron tunneling between two coupled quantum dots (CQDs), a two-state quantum 
system (qubit), using a low-transparency point contact (PC) or tunnel junction as a detector (environment) continuously 
measuring the position of the electron. Here /il and /ir stand for the chemical potentials in the left and right reservoirs 
respectively. 
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Illustration for difFerent behaviors between unconditional and conditional evolutions. The initial qubit state is \a) 



The parameters are: ^ = 1, £" = 0, 



\X\ = D,, and time is in unit of Q, . (a) unconditional, ensemble-averaged 



time evolution of z{t), which exhibits some oscillation and then approaches a zero steady state value, (b) conditional evolution 
of Zc{t). The qubit starts an oscillation, which is then interrupted by a quantum jump (corresponding to a detection of an 
electron passing through the PC barrier in (c)). After the jump, the qubit state is reset to \b) and a new oscillation starts, (c) 
randomly distributed moments of detections, which correspond to the quantum jumps in (b). 
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Transition from quantum jumps to quantum diffusion. Tlie initial qubit state is \a). The parameters are: ( = 1, 
= TV, — Q, and time is in unit of (a)-(d) are the quantum-jump, conditional evolutions of Zc{t) and 



corresponding detection moments with different ITI/lA"! ratios: (a) 1, (b) 2, (c) 3, (d) 5. With increasing ITI/lA"] ratio, jumps 
become more frequent but smaller in amplitude, (e) represents the conditional evolutions of Zc{t) in the quantum diffusive 
limit. The variable ^{t), appearing in the expression of current through PC in quantum-diffusive limit, is a Gaussian white 
noise with zero mean and unit variance. 
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FIG. 4. Illustration of the quantum Zeno effect. Both conditional (in solid line) and unconditional (in dashed line) evolutions 
of the population difference for different ratios of (a)(rti/n)=0.04, (b) 2, (c) 8, are shown. The initial qubit state is |a). The 

other parameters are: (" = 1, £^ = 0, 6 = tt, |T|^ = 20f2, and time is in unit of (2f2)~^. Increasing (Fd/fi) ratio increases the 
period of coherent oscillations between the qubit states, while the time of a transition (switching time) decreases. 
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-6 -4 -2 2 4 6 

FIG. 5. A plot of the noise power spectrum of the current, normahzed by the shot noise level for different ratios of 
(a)(rd/n)=0.04, (b) 2, (c) 8. All the parameters are the same as the corresponding ones in Fig. ^. For small (Fd/O,) ra- 
tio, two sharp peaks appear in the noise power spectrum, as shown in (a). In (b), a double peak structure is still visible, 
indicating that coherent tunneling between the two qubit states still exists. In the classical, incoherent regime Td > 4:Q, only 
one single peak appears, as shown in (c). 
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FIG. 6. Effect of relative phase on the qubit dynamics. The conditional evolutions of Xc{t), ydt), and zdt) with the same 
initial condition (the qubit being in \a)) and parameters {C, = 1, £ — Q, 6 — -k, \T\^ — 4\X\'^ — 4f2), but different relative phases 
are shown: (a)-(c) for 9 = tt, and (d)-(f) for 6 — cos~^(|A:'|/|T|). The relative phase causes quite different evolutions for Xc{t). 
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FIG. 7. Effect of inefficiency on tlie state purity. The quantum-jump, conditional evolution of the purity Pc{t) for different 
inefficiencies, C = !> 0.6, 0.2 (in solid line), and (in dotted line) are plotted in (a) for a initial qubit state being in a pure state 
\a), (b) for a maximally mixed initial state. The other parameters are: £ ^0,e = -K, \T\^ = 4|Afp = 4n. The purity -preserving 
conditional evolution for a pure initial state, and gradual purification for a non-pure initial state for C = 1 ^^re illustrated. 
However, the complete purification of the qubit state can not be achieved for C < 1. 
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